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The gravitational radiation backreaction effects are considered in the Lense-Thirring approxima- 
tion. New methods for parameterizing the orbit and for averaging the instantaneous radiative losses 
are developed. To first order in the spin S of the black hole, both in the absence and in the presence 
of gravitational radiation, a complete description of the test-particle orbit is given. This is achieved 
by two improvements over the existing descriptions. First, by introducing new angle variables with 
a straightforward geometrical meaning. Second, by finding a new parametrization of a generic orbit, 
which assures that the integration over a radial period can be done in an especially simple way, by 
applying the residue theorem. The instantaneous gravitational radiation losses of the system are 
computed using the formulation of Blanchet, Damour and Iyer(1989). All losses are given both in 
terms of the dynamical constants of motion and the properly defined orbital elements a, e, 1 and 
' The radiative losses of the constants characterizing the Lense-Thirring motion, when suitably 

converted, are in agreement with earlier results of Kidder, Will and Wiseman(1993), Ryan(1996) 
and Shibata(1994). In addition, the radiative losses of two slowly changing orbital elements ^0, $0 
' are given in order to complete the characterization of the orbit. 
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I. INTRODUCTION 

> 

\Q , The determination of the signal template of the radiating companion of a black hole is of great importance for the 
forthcoming interferometric gravitational observatories. An inspiralling neutron star or a several solar-mass black hole 
orbiting a massive black hole, as well as a particle constituent contributing to the radiation of the debris trapped by 
the black hole are notable examples of these companions. Either is the case, one wants to describe the evolution of the 
orbit of a test particle in the neighborhood of a black hole under the influence of gravitational radiation backreaction. 
In a highly idealized picture, the orbit of the particle is a Carter geodesic jj] characterized by four constants of the 
^ motion, the energy, the total angular momentum, the rest-mass of the particle and the separation constant. In the 
1 [ spirit of perturbation theory, the radiation backreaction effects may be taken into account by evolving these constants 
■ of the motion, thus picturing the trajectory of the particle by a sequence of geodesic orbits. 

The history of the orbit may be divided in separate epoches. For example, the first thing that happens to a particle 
^ 1 trapped by a nonrotating Schwarzschild black hole on an arbitrary but distant orbit is that it loses eccentricity due 
to radiation losses, and the radius shrinks until the close region of the last stable circular orbit is approached. There 
one has to replace the adiabatic picture with one valid in the strong field regime. 

The description of a generic orbit in the strong field region of a spinning Kerr black hole is difficult (as has 
been pointed out by Thorne due to the special nature of the separation constant which lacks a straightforward 
geometrical interpretation. Recently, progress in treating the separation constant has been reported by Ori || . 

There is a considerable ease in the description of special orbits such as equatorial or quasicircular, and in fact, 
several papers have obtained the signal templates for such special orbits. In order that any of these special orbits 
can be credited as significant contributor to the gravitational signal, one has to make sure that the Carter constant 
evolves to the envisaged special values. 

Several authors have recently published results on radiative losses by systems of spinning masses. Kidder, Will and 
Wiseman Q formulate the problem by using the Lagrangian (proposed by Barker and O'Connell Q) for bodies with 
finite masses and spins. They and Kidder compute the spin-orbit and spin-spin contributions to the momentary 
radiative loss of energy, linear and angular momenta. In the test particle limit their results have subsequently been 
used by Shibata M, in his study of equatorial orbits, and by Ryan ||[|. By computing the backreaction forces, Ryan 
obtains the power and angular momentum losses for circular and generic orbits. 

In our work, we have taken up the evolution of the orbit in the far field region of the black hole. The trajectory 
of the idealized, nonradiating particle is described by the Lense-Thirring picture ||K|[l~l)| . This picture emerges in the 
test particle limit of the scenario described in [Q. The spin S of the black hole is fixed to the z axis. The orbit is 
viewed approximately as an ellipse precessing in a plane which rotates about S. The motion is characterized by three 
constants: the energy E, the magnitude L and spin projection L z of the orbital angular momentum, and by two 
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slowly changing angle parameters \&o an d 3>o- The Lense-Thirring approximation is linear in the spin terms, thus the 
(additive) Schwarzschild effects may be neglected in this treatment. 

In an earlier paper p^j , we have computed the radiation losses from the massquadrupole tensor in the Lense- 
Thirring approximation, by the orbit smoothing method of Landau and Lifshitz However, one may be rightly 
concerned about the degradation of accuracy of this computation, as a price of the relative simplicity introduced by 
the smoothing method. 

In the present paper, we drop the orbit smoothing, and we include the velocity quadrupole tensor Jik in the 
description of the radiative losses. In the test particle limit, this (formally) corresponds to the 1PN approximation 
||. The rest of the paper is organized as follows. In Sec. 2 we describe the orbital motion in the Lense-Thirring 
approximation. In terms of a set of new angular variables (related to Euler rotations), the equations of motion take 
a particularly simple form: all time derivatives depend on the radial variable r alone. 

In Sec. 3, we explore the advantages of having a pure radial equation of motion. We find the turning points at 
f = which allow us to parameterize the radial motion by trigonometric functions. Integrations over the time can be 
replaced by integrations over the parameter, which are straightforward in principle, but can be messy in practice. We 
evade such complications by finding suitable parametrizations. The first parameter £ we introduce is a generalization 
of the eccentric anomaly of the Keplerian orbits. It is employed for computing the half period defined as the time 
spent by the particle between consecutive turning points. For the purpose of averaging the radiative losses, and for 
determining the change in the new angular variables, a generalized true anomaly parameter x is properly defined. In 
terms of this parameter, when integration by the residue theorem is performed, the only pole of the integrands will 
be in the origin. 

In Sec. 4, we compute the rates of change of the constants of motion due to radiation backreaction, to first order 
in S, by employing the radiative multipole tensors of Kidder, Will and Wiseman 0, which originate in the Blanchet- 
Damour-Iyer formalism Jl3| , pT[ . As we are interested in the cumulative effects of the secular motion, we compute 
averaged losses over the period of radial motion. We give the losses in an unambiguous manner in terms of the 
constants of motion E, L and L z . 

As has been pointed out by Ryan 0, the definition of the orbit parameters a and e is subject to ambiguities. In 
Sec. 5, we define the semimajor axis and eccentricity. With advantages explained in the main text, our definitions 
differ both from Ryan's and the ones used in the Lense-Thirring paper. We rewrite all losses in terms of our orbital 
parameters t, a and e. 

In the concluding remarks, we complete the characterization of the changes caused by gravitational radiation by 
giving the change in the ^o^o angle parameters over one orbit. We then compare our results with those of other 
authors. 

We wish to make some comments on the expansion procedure we have adopted. Earlier works employ either the 
expansion parameter e ~ v 2 ~ m/r or the inverse 1/c of the speed of light. Since we eliminate v 2 and f 2 by use of 
the first integrals of the motion, we choose to retain G and c for bookkeeping of the orders in our formalism. The 
correction terms in the Lagrangian, in the equations of motion and radiation losses will carry an extra 1/c 2 factor. In 
our computation, we keep only the terms containing the spin S from among these correction terms, thus we may as 
well choose S as our expansion parameter. 

Our treatment of radiation backreaction effects can be generalized beyond the Lense-Thirring approximation. In a 
follow-up paper |l5|] we will apply our averaging method and parameterization for computing the backreaction effects 
on a bound system of two finite masses. 



II. THE ORBIT IN THE LENSE-THIRRING APPROXIMATION 



We consider the orbit of a small particle of mass \x about a body with mass M and angular momentum vector S. 
In the Lense-Thirring approximation [ ^p[ , the Lagrangian of the system has the form [ flT| 

£= M£! + GmM (21) 
2 r 

where r = |r|, a dot denotes d/dt and the term describing (perturbative) rotation effects is 

E = ^S.(fxr). (2.2) 

The Cartesian coordinates of the test particle are r = {x, y, z}, with origin chosen at the black hole. 
Hence the orbital momentum 
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L = r x p 



(2.3) 



and Runge-Lenz vector 



A= P xL _GMfr 



(2.4) 



satisfy the equations of motion, respectively 



L=^SxL 



(2.5) 



2G 



S x A 



6G 

c 2 r 5 fj, 



(S • L) r x L 



(2.6) 



Here p = dC/dr is the momentum of the orbiting mass. Note that the angular momentum L undergoes a pure 
rotation about the spin with angular velocity: 



(2.7) 



From (2.3) and (2.4), the vector r can be written in terms of the basis {A , L x A} in the plane with normal L: 

L 2 - G[i 2 Mr r • p 



liA 2 



liA 2 



L x A 



(2.8) 



where A = |A|. Inserting this expression of r in (2.6), the direction A/ A of the Runge-Lenz vector is found to rotate 
about a linear combination of the vectors S and L: 



A 



2G_ g _ 6G(L 2 - G^Mr) 
c 2 r 3 c 2 fi 2 A 2 r 5 



A 



(2.9) 



Thus the direction vector rotates about S with the angular velocity £1$ (as does L) and about L with the angular 
velocity: 



6GL(L 2 ~ Gfi 2 Mr) 

^2^5 ( S ' L ) 



(2-10) 



Although the orbital momentum and the Runge-Lenz vector are not conserved (unlike in the zeroth order in S), 
here also there are constants of motion. Because of the fact that the potential is axisymmetric and stationary, the 
orbital momentum component L z and the energy E — rdC/dr — C are conserved. Since the motion of L is a pure 
rotation, its magnitude L is also constant: 



E = L z = L = 



The constants of the motion are 



E = + r 2 (6 2 + sin 2 9 ^)] - 

2 ' r 

L z = fir 2 sin 2 6 cp- sin 2 9 

c 2 r 

L 2 = fi 2 r\9 2 + sin 2 9 ^ 2 ) - 4^^r sin 2 9 £ 



(2.11) 

(2.12) 
(2.13) 
(2.14) 



(Here we are using polar coordinates, x = r sin 9 cosy, y — r sin (9 sin y, z = rcos9). As a consequence, the angle l 
subtended by L and S is also conserved: 



COSi 



L 



(2-15) 
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The first integrals of the equations of motion follow from (2.12)-( 2.14) by simple algebra (we need not recourse to the 
Hamilton- Jacobi formalism). Thus the equation of the radial motion takes the first-order form, decoupled from the 
angular degrees of freedom: 



fi 2 r 2 



GM E GL Z 

! + 2--4^ r -^S , 

r n crfir^ 



the equations for the angular variables 9, ip, however, are coupled among themselves and with the r motion: 

L 2 A cos 2 t\ . L z 2GS 



jj, 2 r A 



1 



sin 



\ir 2 sin 2 i 



(2.16) 



(2.17) 



There is an alternative description of the motion in the Lense-Thirring approximation which bears a more intimate 
geometrical relation to the picture provided by perturbation theory. We shall introduce this formalism by performing 
time-dependent Euler rotations about the origin so as to place the Kepler ellipse in its momentary orientation: 



r = R z (<S>)R x (l)R z (V)v . 



(2.18) 



Here 




(2.19) 



is the initial position of the particle on the plane perpendicular to L (although L is different from Ljy = fir x v, the 
condition r _L L is fulfilled). A rotation about the Cartesian x and z axes is given, respectively, by 



/I 

RxW — C0Si 
\ sin i 




/ cos $ - sin $ " 
R z ($) = sin$ cos$ 
V 1, 



The coordinates of the particle are written in terms of the Euler angles $, I and \I/ as 



x 

V 

z ■ 



r(cos $ cos ^ — cos i sin $ sin 'J') 
r(sin $ cos 'f + cos i cos $ sin 'I') 
r sin l sin ^> . 



Hence we get the relations 



sin l 



cos*vI/ 



(p = $ 



sin 2 6 



sm 



= 1 — sin 2 l sin 2 ^ , 



which allow us to express the time derivatives of the new variables: 

L 



* = 



$ = 



2GS 



(2.20) 



(2.21) 
(2.22) 

(2.23) 
(2.24) 



Note that $ = tts and 1 £ r carries the interpretation of angular velocity of rotation about S and in the plane perpen- 
dicular to L respectively. A n im portan t fea ture of the new coordinates 'J', $ is that their deri vativ es appear linearly 
in the equations of motion ( [j.23| ) and ( [2.24 ) . This is not the case with the polar coordinates ( 2.17 ) because the root 
of 9 cannot be extracted. The geometrical interpretation of the angle variable \& assures that it is monotonously 
changing. We choose '5 to increase with the motion evolving. 

There is a freedom in the precise way one introduces the three Euler-angle variables in place of the polar angles. 
Even though the equations of motion are made simple by use of the angular variables ( 2.19 - 2.21 ), the description of 
the momentary plane of orbit is tied to a different choice of the Euler angles (^', J, where the initial z axis is 
chosen along L^v rather than along L. These angular variables appear (with the appropriate change in notation) in 
the Lense-Thirring paper. The angle if is not a constant, and the resulting equations are less simple. The primed 
angles will not be employed in our computations, and we relegate their discussion to the Appendix. 
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III. PARAMETRIZATION OF THE ORBIT 



We now proceed with the solution of the equations of motion. First consider the radial equation ( [2.16 ). For the 
turning points r — 0, we get the cubic equation 



2Eur 3 + 2GMu 2 r 2 - L 2 r - A%uL z S = 

c l 

In the no-spin limit, S — 0, the solutions are (the unphysical) r = and 

-GMu =F A 



r± 



2E 



where Ao is the length of the Runge-Lenz vector to zeroth order 

A = Jg 2 M 2 u 2 + 



2EL 2 



li 



(3.1) 



(3.2) 



(3.3) 



When a small spin term is present, the roots of (3.1) will have a slightly different form from (3.2): + e and r± + e±, 
where e and e± are assumed to be small. Note that the constants E, L and L z will not coincide with their zeroth 
order limits. The constant A is still defined by (3.3), although it is not identical with the length A of the Runge-Lenz 
vector (which is not constant). Writing the cubic equation as 2E/i(r — r + — e + )(r — r_ — e_)(r — e) = 0, we get 

2En[r 2 - r(r + + r_) + r + r_]r + 2Efi[-r 2 (e + + e_ + e) + r(e + r_ + e_r + + e(r + + r_)) - er + r_] = . (3.4) 

The first term is the zeroth-order part of the original equation, written in terms of the roots, thus it equals 
(AG/c 2 )/iL z S. Inserting the chosen form of the roots we get 



4G^L Z S 
c 2 L 2 ' 



2G[iL z S 
c 2 L 2 



(^" 



(3.5) 



The orbit is supposed to lie at large distance, thus the first root will not occur. 

Using the above results, a parameter £ generalizing the eccentric anomaly of the Kepler orbits is readily introduced: 



G/iM 2G/iL z S 



2E c 2 L 2 

such that r— is at £ = 7r and £ = 0, respectively. From here: 

/ A 2G 2 pi 2 ML z S 



A Q 2G 2 n 2 ML z S 
2E + 



dr 
rf£ 



2E 



c 2 A L 2 



- 2 A L 2 



sin£ 



cos£ 



(3.6) 



(3.7) 



For equatorial orbits, this parametrization has been employed by Shibata J7j. Expressing 1/f from (2.16) and lin- 
earizing in S, the expression for dt/d^ is found: 



dt 



d£ fd£ 



1 dr n 2 (GMn - A a cos£) 2G 2 M^L z S cos£ 



(-2/^)1 



(-2fiE)iL 2 A c 2 



Integration of (3.8) from to 2tt gives the orbital period: 

GM/j? 



T = 2ir- 



(-2/*S)» 



(3.8) 



(3.9) 



However the parametrization ( |3.6| ) is inconvenient for carrying out the averaging of t he l osses, because it yields a 
complicated array of poles of the integrands. Hence in the sequel we will, in place of ( |3.6| ), find a parametrization 
r = r(x) satisfying the following two criteria: 



(a) r(0) = r„ 

(b) dT 
d(cosx) 



and r(ir) = r max 
-(7o + Sji)r 2 , 



(3.10) 
(3.11) 
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where 70, 71 are constants. That is, we want to keep the properties of the true anomaly parametrization of the Kepler 
orbit. Property (b) generalizes Kepler's second law for the area fllf l . The unique parametrization satisfying both (a) 
and (b) is: 



L 2 



4GL Z S A (2G 2 M 2 ^ + EL 2 ) + GMfi(2G 2 M 2 ^ + 3EL 2 ) cosx 



fx(GM(j, + A cosx) A L 2 c 2 



(GMfj, + A cosx) 2 



Hence 



Mo 4G 2 M^L z S.„ n/f2 



V- 



and to first order in S: 



dr 
dx 

dt 1 dr fir 2 
dx r dx L 



A L & c 2 



(2G^Af V + 3EL 2 ) 



r smx 



1 



cosx) 



We will need also the relation of the type \& = ^(x)- This can be obtained from integration of ( p.23| ) 

2Gfi 2 L z S 



* = * 







X ■ 



2 L 4 



-(3GMnx + Ao sm X ) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



The constant of integration ^f measures the angle subtended by the semiminor axis (x = 0) and the node line 
at \& = 0. In a perturbative picture, the orientation of the momentary ellipse of the orbit is reinterpreted at each 
revolution, and the valu e of t he integration constant is subject to a corresponding shift. 

Our parametrization (3.12) makes especially simple the integration over one period of all expressions F = F{r, ty) 
conta ining in the d enomi nato r no \f r dependence and no other r dependence than r 2+n , with n a positive integer. Note 
that (2.1C), ( [2.23| ) and (2.24) have this property. Further, as will be shown in the next section, the instantaneous 
radiative losses of the constants of motion share this property. In all of these cases one has: 



F(r(t),^>(t))dt 



%(r(x),*(x))^X, 
dx 



where the r dependence of the denominator of the integrand is especially simple: 



< F >- 



F(r(t),*(t))dt 



(3.16) 

The average over one period is 

(3.17) 



with T given by Eq. (|3J). 

The integrals may conveniently be evaluated by use of the residue theorem. One introduces the complex variable 
C = e* x , and the integral is taken over the unit circle in the £ plane. The value of the integral is given by the sum 
of the residues of the poles inside the circle. In all relevant cases our parametrization assures that the only pole is at 

C - 0. 



Simple examples are the integrations of (2.24), (2.23) and (2.10): 

,2G 2 Mfi 3 



2irS- 



= 27r - 2ttS 
An L = - 2nS 



c 2 L 3 

6G 2 Mn 3 L z 
6G 2 M/i 3 L z 



2 L 4 



(3.18) 
(3.19) 
(3.20) 



Here ( 3.18D is the precession angle of the orbital angular momentum L about S. The interpretation of the variab le \& 
as the polar angle in the plane perpendicular to L is in accordance with the identical first order terms in ( 3 . 1 9|) an d 
(3.20). The ellipse is found to counterrotate (for L z > 0). The periastron shift is a combination of (|3.1 <j) and (3.20). 
Both and A<i> will be unchanged in the averaged-motion approximation of Landau and Lifsitz |11|J12[| . 

Our angles are related to the polar angle (p for equatorial orbits (1 = 0) by ip = $ + ^. This relation helps 
comparison with Shibata's expression (2.19) for the periastron shift Aip. 

Ryan || has introduced a parametrization r = r[ipRj which, when rewritten in terms of the constants of motion 
E, L and L z (employing his definitions for the orbital parameters), will read: 



G 



AGL Z S A a (2G 2 M 2 fi 3 + EL 2 ) + GMu(2G 2 M 2 fi 3 + 3EL 2 )cosip R - (Mo/2) sin 2 tjj R 



fi(GMu + Aq cos ip R ) A a L 2 c 2 (GMu + A cos ip R ) 2 

(3.21) 



This differs from our parametrization (3.12) by the last si n 2 ip R ter m, wh ich shows that Ryan's parameter ip R and 
our x reach the turning points simultaneously. In fact both ( |3.6| ) and ( 3.12| ) can be supplemented by arbitrary terms 



involving a sin squared factor, and they still will satisfy (a). However none of these parameters, and among those 
Ryan's, will satisfy the generalized second Kepler law (b). 

IV. INSTANTANEOUS AND AVERAGED LOSSES 

The mass and current quadrupole tensors of a test particle with coordinates r,; = (x,y, z) are 

Iii = n(r i r j ) 3TF (4.1) 

J l3 = -„[r,(r x §)/ TF + ^{nS^ , (4.2) 

where STF means symmetrization and trace removal in the free indices. The rates of the energy and the angular 
momentum losses in the test particle 1PN approximation are, respectively 

dE _ G / d 3 Ijj d 3 Ijj 16 ePJjjdPJjj } 

dt 5c 5 V dt 3 dt 3 Abe 2 dt 3 dt 3 ) 1 ' ' 



dLi G / 2 d 2 1 j r d 3 Ij~ r 32 d 2 J j r d 3 J) 



/ z ™ i — 1 "3 r 1 " kr ) (A 41 

dt ' c 5 ljk \5 dt 2 dt 3 Abe 2 dt 2 dt 3 1 ' 



Evaluating the right-hand sides by use of (2.16), (2.23) and (2.24) and dropping higher-order terms in S, the total 
power is 

dE = _ , 2 2 a) 8 SG 3 ML cos l , ^ _ a 2J a . 

df 15 c 5 r 6 v ; 15 /ic 7 r 8 v ; y ' 

This is independent of the angles \& and The loss in L is evaluated by 2L5L = 5L 2 = S(LiLi) = 2LiSLi and it 
has the form 

dL =+ 8G2 LM , 2 E 2 _ 3 j2\ + S G 2 S cos i GEM 2 r 2 + 3 Q 2 ^M 2 r - 11 MGL 2 [i + 18 EL 2 r) . (4.6) 
dt OjUcV y ' 15 c'fj,r° 

The momentary loss in L z is 

^=+ 8G2LM ™ SL (2Eur 2 -3L 2 ) 
A G 2 S 

[( 6 G 2 M 2 P V + 36 n Er 2 L 2 + 6 GL 2 M[i 2 r + 18 L 4 + 24 fi 3 Er 3 GM) sin 2 i/> 

15 fjre'r' 

-39 GL 2 M [i 2 r + 36 L 4 - 36 // £r 2 L 2 ] sin 2 i + 3 rr L/i sin 2^ sin 2 i (2 GAf/i 2 r — 3 L 2 + 6 En r 2 ) 

-6 G 2 M 2 fi 4 r 2 + 22 GL 2 M[i 2 r - 24 fi 3 Er 3 GM - 36 // £r 2 L 2 } . (4.7) 



After parameterizing by % and averaging over the period T as described in the previous section, we obtain 

dE ^ - 4 (rEV?' 2 &M £2lA + EG 2 L 2 M 2^ + 425 g^^h 



eft / 15 c 5 L 7 V2 



2 S cos lG 2 { y 3/2 (520 ff 3£ 6 + 1074Q G 2 E 2 L i M 2 3 + 24990 G i EL 2 M i u 6 + 12579 G 6 AfV) 
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dL\ 16 G 2 {-Eaf 2 M , , . 



eft / 5c 5 L 4 V2 



-(UEL Z + 15 G M n ) 



4 G 2 (-En) 3/2 S cost 



(1188 £ 2 L 4 + 6756 G 2 EL 2 M 2 V ? + 5345 G 4 M 4 fi 6 ) 



l5c 7 L 7 V2 

1\/T pns / 

(14££ 2 + 15 G 2 MV 3 ) 



\ <ft / 5c 5 L 4 V2 

4G2 ^(-^) 3/2 (n88£; 2 L 4 6756G 2 SL 2 M 2 3 + 5345 G 4 M 4 6) (4 g) 

15c 7 L 7 \/2 

- 2G2 £(-^) 3 (_ 2 sin 2 t (n72 ^4 + 5892 G 2 EL 2 M 2^ + 4325 G 4 M 4 m6) 

5 c 7 L 7 v 2 

8 G 2 ^-^ cos(2* ) sin 2 1 4 ^2^2 3 gg^.j 

5c 7 L 7 V2 

Among these losses averaged over one period of the radial motion, the only quantity depending on the initial angle 
^0 is {dLz/dt}. This is to be interpret ed s uch that the rate of loss of the angular momentum component L z varies 



with the position of the periastron. In (4.8), the terms with cos^^o) average to zero when the precession time scale 



is short compared to the radiation reaction time scale. (For a more extensive discussion of this subject, cf. ]9j].) 

V. AVERAGED LOSSES IN TERMS OF ORBIT PARAMETERS 

The major axis a and the eccentricity e of the Kepler motion are constants and together with 1 they determine the 
full set of the constants of motion E, L and L z . In the Lense-Thirring picture of a perturbed Kepler motion E, L and 
L z contain first-order terms, and the corresponding constants a and e should also contain first-order terms. There 
is an ambiguity in how to introduce these orbit parameters, as noted by Ryan He chooses these parameters by 
referring to the quasicircular orbits in the Kerr metric. This method is passing when wishing to approximate Carter 
orbits. However, there appears to be no reason for invoking phenomena tied to black holes when describing arbitrary 
weakly bound spinning bodies. In fact, we want to keep our treatment general enough to cover all axially symmetric 
gravitational fields with first order contributions from the spin-orbit interaction term in the Lagrangian ( |2.2| ). Thus 
we choose as the definition of the semimajor axis a and eccentricity e: 



which implies: 



a(l±e), (5.1) 



„_ GM,l 2GSCQSC 

2a ^ c 2 a^GM(l-e 2 ) J 

r2 2 /, 2n ( ■> 2G5cosi 3 + e 2 \ , . 

L = GMfx a(l - e 2 ) 1 - ___=__- . (5.3) 



c 



2 ai v /GM(l-e 2 ) 



1 -e 



These relations can be inverted: 



GMfi f 4EL Z S\ 
2E V ~ c 2 ML 2 J 



1-e 



2 



2EL 2 



G 2 M 2 n 3 

The averaged losses in terms of a and e then read: 

/ dE\ , , x G 4 M 3 u 2 

- (37e 4 + 292e 2 + 96) : 



1 + S (G2MV + ^ 2) 



(5.4) 
(5.5) 



dt / y 'lb(l-e 2 y/ 2 a 5 c 5 

+ (491 e 6 + 5694 e 4 + 6584 e 2 + 1168) 1" 7 ' / 2 ' (5.6) 

30 c'a L6 i z {l — e A y 



8 



dL 

dL z ' 
~~dT 



-(7e 

- (7e 2 + 8) 



4G 7/2 M 5/ 2/i 2 

5a 7 / 2 (l-e 2 ) 2 c 5 



(549 e 4 + 1428 e 2 + 



G 4 SM 2 fi 2 cost 
15c 7 a 5 (l-e 2 ) 7 / 2 



(5.7) 



, ^G^M^V cos/ 



5c 5 a 7/2( 1 _ e 2)2 



(285 e 4 + 1512 e 2 



G 4 SM 2 fi 



2, ,2 



30c 7 a 5 (l-e 2 ) 7 / 2 



+ (461 e 4 + 1456 e 2 



G 4 SM V cos 2 1 
10c 7 a 5 (l-e 2 ) 7 / 2 



(13 e 2 + 20)2e : 



,G A SM 2 fi 2 sin 2 tcos(2^ ) 
5c 7 a 5 (l-e 2 ) 7 / 2 



Differentiating (2.15) one has: 



di 
dl 



1 



L sin t 



cos L{ 



dL 
~dt 



dL z 
~dt 



(5.8) 



(5.9) 



Those terms in the brackets which are zero-order in S. cancel. 



The averaged losses of the ellipse parameters can be obtained by differentiating (5.4) and (5.5) and from (5.9) 

G 7/2 M 3/2 eM 



de 
di 

da 
~dt 



G M fie 



Yl\e A + 304 



di 
dl 



— ) = S sin 



15 (1 - e 2 ) 5 / 2 a 4 c 5 

64G 3 MV /37 4 73 
5o 3 c 5 (l-e 2 ) 7 / 2 ^ 96 6 + 24 f 

G 7/2 M 3/2 M 



(1 - e 2 ) 4 a n / 2 c 7 



/ 1313 4 932 2 
Scob^— e + — e 



64 G 7 / 2 M 3 / 2 /iS cos l ( 121 6 j 585 
T a 9 / 2 c 7 (l - e 2 ) 5 l"64~ e ' 



32 



1172 \ 

124 , 133 S 

e 2 H 

3 12 



l ll/2 c 7( 1 _ e 2)4 



,„ T J 26 4 „ 2 \ 19 4 252 2 244 
cos(2* o) -e 4 + 8e 2 +ye 4 + — e 2 + — 



(5.10) 
(5.11) 
(5.12) 



Let us take an example of a binary system contaning a massive black hole with mass 100M Q and spin S » 0.6M 2 , and 
where the companion's orbit has the semimajor axis 10 4 km, i.e., one-hundredth of the Hulse- Taylor system. Then the 
increase in the angle l per loss of logo (a dimensionless quantity) is < di/dloga >= (G/m/ a 3 ) 1 ^ 2 S sin l 10 -3 sint 
for a circular orbit. 

The enhancement factor for an eccentric orbit is 



/(e) 



285e 4 + 1512e 2 + 488 - 6(26e 4 + 40e 2 ) cos(2W ) 
4(37e 4 + 292e 2 + 96)(1 - e 2 ) 1 / 2 



(5.13) 




FIG. 1. The enhancement factor. 
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The enhancement factor is monotonously increasing (Fig.l) from /(0) = 122/121 to make the change in i comparable 
to that in log a only when the eccentricity approaches the unit value to one part in 10 6 . 



VI. CONCLUDING REMARKS 



The Kepler orbits are completely characterized by the energy E, the angular momentum vector L and the Runge- 
Lenz vector A at the initial time tQ. These are six independent constants of the motion. In the Lense-Thirring picture, 
the constants of the motion alone do not suffice for the full description of the orbit. We complete the characterization 
by specifying the Euler angles "Jo and $o of the periastron . Th eir evo lution over one period, under the influence of 
the spin (disregarding radiation backreactio n) is given in ( 3.1 8| ) and (3.19). Taking into account the gravitational 



radiation backreaction, the shifts ( |3.19 ) and ( |3 .18 ) will vary from one period to another, as the 'constants' L and L 



evolve. One may interpret these as the radiation- induced changes in the parameters "to and $o: 

i n,r . .3 

5^o ■= A*(T)-A*(0) = 2nS 2r / SL Z (6.1) 

C Li 

(5$ := A$(T) - A$(0) = -2nS — jy^-SL , (6.2) 
where SL = TidL/dz) and we have used tha t in z eroth ord er SL Z = 5L cos t. The total shift in the Euler angle of 



the periastron point is given by the sum of (3.19) and J6.1| ), and likewise for <i>o. 

We continue with some comments on the radiative losses in the constants of motion. The terms independent of 
S agree with the results of Peters and Mathews P^pl ■ The terms linear in the spin S in our averaged losses agree 
with Ryan's corresponding results, when the parameters a, e and t are converted to his parameters. In the present 
paper, the energy and angular momentum losses have been directly computed. This is to be contrasted with Ryan's 
method where the backreaction forces have been used. The perfect agreement achieved underscores the reliability of 
these computations. 

When l = 0, the averaged losses in L are equal to the averaged losses in L z up to first order in S. Hence the 
equatorial plane of the orbits is stable with respect to radiation losses. 

As we have shown in the previous section, the inclination of the orbit cannot change significantly in the adiabatic 
era. Hence we conclude that it is insufficient to restrict the computation of signal templates from a black hole - test 
particle system to equatorial Carter orbits. 

In the limiting case of circular orbits, e — > 0, there is no distinguished point on the ellipse. The dependence of 
(dL z /dt) on Vl'o vanishes. As expected from the axial symmetry, the radiative losses, for arbitrary excentricity, do 
not depend on the initial angle $o- 

In the present paper, we have developed a toolchest for the computation of radiation losses, which include a new 
parametrization of the orbit, new angular variables and the application of the residue theorem for obtaining time 
averages. Currently we are employing these procedures with the inclusion of finite mass effects. |15j 
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APPENDIX A: ALTERNATIVE ANGULAR VARIABLES 

Introduce three angle variables tl and where *S?' is the angle of rotation about Ljv, bringing the position of 
the particle to the node; i! is the angle subtended by Ljy and S and $' is the angle of precession of Ljv about S. The 
angle is the argument of the latitude, i' the inclination and $' the longitude of the node |lC| ] . Unlike t, the angle i 1 
is not constant. Let us express if in terms of t and Computing the z component of Ljy in two different ways we 
get the two sides of the equation 

{ /i 2 r 4 (6 2 + sin 2 9 <p 2 ) }2 cosi' = fir 2 sin 2 6 >p . (Al) 



With the time derivatives of polar angles (2.17) after linearization in S, the previous equation gives 
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COS l' = COS i(l + 



2GnS 

c 2 L z r 



sin tcos^ *) . 



(A2) 



Here we are using that the primed angles are equal, to zeroth order in S, with the corresponding unprimed ones. 
Taking the time derivative of i! we get 

2GS sin i cos * 



, , - (/irf cos \& + 2Lsmty) . 
c Lt 

With the new time-dependent Euler rotations the coordinates of the particle become 

x = r (cos $' cos ty' — cos i' sin <&' sin 
y = r(sin <&' cos + cos tf cos <E>' sin 
z = r sin i! sin 

and the time derivatives of the polar angles are 

• sin t' cos ^' •, cosi'sin^E'' ., „ .■>,.•> -r i 

9 = — — l , sir 9 = 1- sur i sin 2 V 

sin sin 



(A3) 



(A4) 



(A5) 



• , cos i' ■ . sin J sin 2'$)' , 



sin 2 6> 



2 sin 2 6 



(A6) 



Inserting these equations in (2.17), we get the changes of the primed Eulcr angles in terms of the unprimed ones 



• , L GS cos i .„ r » 

* = — - H „ (2L - fxrr sin 2*) 



2 Lr 3 



.. 2G5sin* , . T or . T , 
$ = 5 — 5 — i\irr cos * + 2L sin *) . 

Integrating these equations, we get the new angle variables and their shifts over one period 

. » AT/ „ „ ry 6G 2 M/i 3 L z A . „ 2G 2 M/i 3 
At' = 0, A*' = 2tt - 2^5 ^ , A$' = 2ttS — j-j- . 

C J-j C -Lj 

Note that these changes in the primed variables are equal to the shifts (3.1£) and (|3.1S|) of the unprimed variables. 



(A7) 
(A8) 

(A9) 
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